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THE  SUPERSONIC  FLOW  AROUND  SYMMETRICAL  TORS I OWED  AND 

CURVED  DELTA  VINOS,  TAKING  INTO  CONSIDERATION  THE 
SEPARATION  OF  FLOW  AT  THE  LEADING  EDGES 


by  ELIE  CARAFOLI  and  STEFAN  STAICU 

Institute  of  FLUID  MECHANICS 

Academy  of  the  Socialist  Republic  of  Romania 

In  this  work  the  supersonic  flow  around  deformed  thin  delta 
wings  1 8 studied,  having  the  distribution  of  Incidences 
symmetrical  and  varying  proportionally  with  Xj  and  x2. 
Taking  into  consideration  the  separation  of  flow  at  the 
subsonic  leading  edges,  the  distribution  of  pressure  and 
and  aerodynamic  characteristics  of  the  wing  are  determined, 
through  the  intermediary  of  an  imaginary  wing,  equivalent 
from  an  aerodynamic  point  of  view,  with  the  real  wing* 

1.  PRELIMINARY  CONSIDERATIONS 

In  the  present  work  we  will  study  the  supersonic  flow  with 
the  separation  of  flow  at  the  leading  edges  of  deformed  thin 
delta  wings,  in  such  a way  that  the  incidence  is  symmetrical  in 
reference  with  the  axis  of  symmetry  and  varies  proportionally 
with  Xj  and  x2  (fig*  1 ) • 


As  with  the  plane  thin  delta  wing,  with  constant  incidence 


1 
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(1),  (9)*  in  this  case  the  flow  separates  at  the  leading  edges, 
creating  as  well  a vortex  layer  which  winds  itself  approximately 
in  the  fora  of  two  horns,  which  are  then  transforaed  for  the 
no8t  part  of  the  vortex  generating  intensities  into  two  concentrated 
vortex  nuclei,  situated  syaaetrically  in  reference  with  axis 
Oxj  (fig.  1)  and  defined  by  the  coordinates  c and  t.  The  system 
of  two  concentrated  vortexes  considered  will  bring  changes  on 
the  flow. 

The  resulting  notion,  which  becomes  more  complicated  now, 
will  be  studied  on  the  basis  of  the  theory  of  motion  of  the 
second  order.  Towards  this  goal, we  will  follow  the  road  used 
in  previous  papers  (1),  (3)  and  (9),  where  solutions  were  given 
for  thin  delta  wings  with  constant  incidence  and  with  anti symmetrical 
constant  incidence  respectively  (forced  antisymmetry)  in  reference 
with  0xlt  which  led  us  to  conical  motion  as  a natter  of  fact. 

In  this  manner,  we  will  allow  that  the  effect  of  the  falling 
off  of  flow  at  the  leading  edges  of  the  wing  and  therefore  of  the 
form  of  the  vortex  nuclei  on  the  higher  side  of  the  wing  is  that 
of  changing  the  vertical  velocities  and  axlses  of  disturbance, 
becoming  finite  at  the  edges.  We  can  equate  this  complex  effect 
with  that  of  a wing  with  a variation  corresponding  to  the  incidence 
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or  with  a variation  of  the  vertical  velocities  of  disturbance. 
Through  this  we  substituted  the  horizontal  velocities  of 
disturbance\«eewfr  equivalent  distribution  of  vertical  velocities. 

We  will  obtain  a symmetrical  distribution  of  vertical 
velocities  od  the  higher  side  of  the  wing,  which  will  correspond 
in  this  way  to  an  imaginary  thin  delta  wing,  with  variable 
incidence,  having  at  the  same  time  a finite  velocity  at  the 
leading  edges.  It  will  be  allowed  that  the  real  thin  wing,  which 
has  in  a certain  way  finite  velocities  at  the  edges  through  the 
effect  of  the  separation  of  flow,  is  equivalent,  from  an  aero, 
dynamic  point  of  view,  with  an  imaginary  wing,  having  the  same 
variation  of  incidences  which  we  defined  above. 

Since  the  incidence  is  variable  on  the  surface  of  the  wing, 
the  axis  on  which  the  horn  winds  will  be  a curve,  and  the  vortex 


L 
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generating  Intensity  of  the  nucleus  Is  variable  along  the  axis* 
proportional  with  the  square  of  the  span  of  the  wing.  For 
slmpllf icatlon,  in  the  following,  the  axis  on  which  the  vortex 
nucleus  is  situated  is  considered  a straight  line* 


Noting  further 

V*  — — vm,  u>;  = — ai  Um  (wt  ==  — «a  Vm,  Wj  = — «,  VJ  (I 

X>n  the  higher  side. 

the  vertical  velocities  and  the  incidences  w„,(Xu'^reSperetively# 

U 

/ / 

wi.CXi  on  the  lower  side,  for  the  imaginary  thin  wing  we  can  write 


the  relation 


(*  w,  dy  ■+■  ( wi  dy  = 2 wl  ( v> * Um) 

J»  -o 


where  the  velocity  w on  the  real  wing,  in  conical  motion  of  the 
second  order,  is  the  homogenous  function  of  the  first  order, 

k - ± woiJ*  = ± u>oi y)  = — *i(«io  T HivWm  (*  “ f1)'  (3) 

in  which  the  term  WjqXj  corresponeiSto  the  Incidence  with  natural 

symmetry  (the  curved  delta  wing),  and  a woix2  corr«8Ponds  to  the 
incidence  with  forced  symmetry  (the  torsloned  delta  wing),  the 
sign  (♦)  being  considered  for  the  right  part  and  (-)  for  the  left. 

In  the  same  way  as  in  (1),  (3)  and  (9)«  we  will  split  the 
Imaginary  wing  into  three  wing  components  corresponding  to  the 
distribution  of  vertical  velocities  above. 

1)  The  thin  wing  with  the  variation  of  chosen  corresponding 


Incidence*  In  such  & way  as  to  respect  somewhat  the  significance 
of  the  phenomena  of  pressure  changes  at  the  leading  edges*  In 
this  way  a thin  Imaginary  wing  with  finite  velocities  at  the 
leading  edges*  and  equal  and  opposed  direction  on  the  two  sides* 
higher  and  lower*  Is  obtained* 

2)  The  wing  of  symmetrical  thickness*  having  the  slope 
variable  In  the  same  way  as  the  Incidence  of  the  first  wing.  This 
wing  combined  with  that  from  1)  has  different  pressures  on  the 
two  sides*  as  happens  In  reality. 

3)  The  third  wing  will  have  a symmetrical  thickness*  with 
the  slope  also  symmetrical*  however  In  such  a way  that*  combined 
with  the  wing  from  2l_*  an  average  nought  slope  will  be  obtained, 
corresponding  to  a real  thin  wing. 

2.  DETERMINATION  OF  THE  AXIS  OF  DISTURBANCE  VELOCITIES 

In  continuation  we  will  determine,  for  the  three  Imaginary 
wing  components*  the  axis  of  disturbance  vfe&ocltles*  being 
necessary  for  the  calculation  of  the  distribution  of  pressures 
and  aerodynamic  characteristics  of  the  resulting  Imaginary  wings, 
which  are  presupposed  to  be  the  same  as  those  of  the  real  thin 
delta  wings*  having  the  Incidence  defined  by  (3)* 
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However,  the  motion  around  the  wing  being  conical  of  the 


second  order,  we  will  utilize  the  same  methods  used,  considering 

IJorj 

in  this  sense  the  section  obtained  in  normal  plane  Oyz  imr  the 
direction  of  undisturbed  flow  U^and  having  the  coordinates 


y =~,  z = — , 

*1 


(4) 


axes  Oy  and  Oz  being  parallel  with  0x2  and  Ox^  respectively.  We 
will  further  make  a similar  transformation  with  that  given  b> 
Busemann  (fig.  2): 


1 — JJ«2a  ' 4 1 — (3  — H + ij), 


(5) 


obtaining  a plane  which  has  the  property  of  keeping  the  track 
of  the  wing  (T»y,  z * 0)  in  the  true  magnitude. 


As  we  know,  in  this  plane  the  first  derivatives  of  the 
velocities  of  disturbance  u,  v,  w are  harmonic  functions  and  the 
conjugated  functions  can  be  associated  respectively  in  such  a 
way  as  to  obtain  analytic  functions  of  complex  variables 

* = V + ij.  (8 

We  will  study  further  each  wing  defined  above  in  turn. 


1).  The  thin  wing.  As  a result  of  the  effect  of  the  two 
vortex  nuclei,  the  vertical  velocity  on  the  reml  wing  is  modified, 
as  well  as  on  the  first  wing  component  defined  above.  Thus,  for 
the  points  contained  between  ( -s<y<s  ) from  the  track  of  the 
wing  from  the  plane  x * Y + 1^(6),  the  vertical  velocity  is  the 
same  as  the  form  given  by  (3)* 


— = u>\%'  ± <'  y = — (*io’  =F  < y)  l/.,  (7 

jjj 

where  parameter*  w^q  and  Wqj  correspond  to  the  abscissa  y * s, 
and  for  the  Intervals  (-2,- s)  and  (s,t)  we  will  write 

— = «lo(y)  ± “’ii(y)  y =• -*(«!•  (y)  T ««(y)y)  (*■ 

®i 

so  that  at  the  edges  of  the  wing  we  will  obtain 

— = w[\}±icWl  (* 


The  continuing  variation  of  the  vertical  velocities  ( or, 

more  precisely,  of  the  parameters  w10(y)  and  wm(y))  or  the 

[ continuing/ 

Incidences  corresponds  to  th'e^d 1 s tr 1 but 1 on  of  elementary  edges 
situated  on  the  wing  in  the  interior  of  the  interval  considered, 
which  give  in  each  point  y ■ n.  the  elementary  drop 


(» 


Keeping  in  mind  previous  papers  (2),  (4),  the  contribution 
in  the  expression  of  the  axes  of  disturbance  velocities,  of 
elementary  edges  situated  in  points  y *n  on  the  thin  wing  with 
subsonic  edges  (/<l/B),  through  the  application  of  similar 
hydrodynamic  methods,  will  be 


W = (q(£  + *?«’)  c°8h-‘ 


l/d  + *))(!  — *) 
\ 2i  (T)  — x) 


d >)• 


(11 


In  a similar  way,  for  the  contribution  of  elementary  edges 
situated  on  the  left  of  the  origin,  y ■ -•'i,  we  will  have 


d«l"  = (ftf  + a-  qti)  co8h->  1 d* 
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These  contributions  of  the  edges  In  the  expressions  of 
axes  of  disturbance  and  vertical  velocities  on  the  wing  are 
achieved  placing  some  singularities  of  the  second  order  (In 
this  case  - sources)  on  the  track  of  the  wing  In  plane  x*Y+i^(6). 

In  this  way  we  succeed  In  acquiring  In  the  calculations 

the  effect  of  the  vortex  sheets  which  fall  off  from  the 

edges  and  of  the  considered  concentrated  vortex  nuclei,  through 
icontlnulngj 

out  the'dl s tr lbut  1 on  of  sources. 


Next  we  will  allow  a simpler  form  of  division  of  sources, 
which  will  satisfy  conditions  Imposed  by  the  problem,  for  obtaining 
concommltant  axes  of  disturbance  and  vertical  velocities  Indicated 
above  on  the  basis  of  observations  and  conforming  with  experimental 
results. 

Thus,  we  will  chose  a llnlar  variation  of  intensities  of 
the  sources 

9'  - («““»  + 9,(/')  (i  — y) . (13) 

where  the  constants  q*^n^and  q*0)  correspond  to  natural  symmetries 
forced  respectively,  of  the  incidences.  Since  the  flow  Is  conical 
of  the  second  order,  as  seen  in  (11),  (12)  we  will  write 

qm—  “"j»  = ^1 (14) 
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keeping  In  Bind  these  divisions  of  sources,  starting  from 
(11)  and  (12),  we  will  consider  the  contributions  of  all  the 
elementary  distributed  edges,  as  well  as  t<n  the  subsonic  edges 
(2),  In  order  to  obtain  the  following  expression  of  the  axes  of 
disturbance  vleocltles  of  the  first  wing  components: 

- Tif=r-' + IS  to“+9“‘|c“b"  + 


+ *^- (?i>— • «ii«)cosh  ']j 


jl  ~t~  T') ) {l  ~t  "•*>) 
2 i ( r)  + x) 


Introducing  q2o  and  421  ln  (I**)  and  making  the  calculations 


we  obtain 


A<,)  -4f),OHh',!(^i  ■* 

+ «*  - x)  (« + •)  (1  - ^r) t0Hlrl  ^ + 

+ [«So  *({/  1—^  — 2 COB-*  -yj  + gi,  X*  COH  1 -* j j/l  — , (16) 


ln  which  Ago*  A22*  q20»  421  are  8ome  constants  which  will  be 


determined  below. 


2).  The  wing  of  symmetrical  thickness,  with  the  slope  equal 
with  the  Incidence  of  the  first  thin  wing.  Through  the  Introduction 
of  this  wing  of  symmetrical  thickness  the  accentuated  peaks  of 
pressure  on  the  lower  side  of  the  wing  are  removed,  where  the 
distribution  of  pressure,  obtained  through  the  superposltlonlng 
with  the  thin  wing  1).,  Is  different  from  that  on  the  higher 
side.  Following  the  general  method  of  conical  motion  (2),  (4), 
corresponding  to  a wing  of  symmetrical  thickness,  with  the  slope 
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defined  by  the  same  distributions  of  sources  given  by  the 
relations  (13)  and  (14),  we  will  write  for  the  axis  of  disturbance 
velocity  the  following  expression: 


CUXt  = ~ [ (<ho  + 9*.  X)  COHh-‘  ]/ 


— Bx) 
Tx) 


drj  -f 


*■ + L’  m 

which,  in  the  course  of  the  accomplishment  of  the  calculations 
becomes 

i i / «. 

— BHx 


— = -L  |(9*u  + q*n  ®)  £(l  — a;)  |i  _ ~~j  cosh-1  - 

_ (*  _ ,)  (i  _ i±JE)  C08h-i  , 

l 21  ) B{»  — a:)]  ^ 

+ (9»o  — 9?i  |(i  + x)  |l  — I* 2i~  j cos^~l  ~ 


B(l  — a:) 


+ B'lx 


.(8  + JC)(1„lZl£)cotlh-lL±J«li?] 
V 21  ) B(»  -f  x)J 


21  ) b(l  + x) 

+ B*»x  1 


+ 


+ 


B(*  + x) 

~t  [2/  (ain-1  Bl  — sin ~ 1 Hs  + £L~  B'  l\~  Kl  — ^ «*)  .. 


2BI 


— (sin-1  Bl  — sin-i  ha)  x ? /1  — B*  x*j  + l, 


-)*- 


(18) 


where  the  term  L is  due  to  the  slope  of  the  leading  edge: 

v<i> 

— “Yft1  + Ytf*.  (W 

3).  The  wing  of  symmetrical  thickness  compensating  for  slope. 
Through  the  effect  of  the  wing  from  point  2).,  the  resulting 
wing  became  the  wing  of  "symmetrical  thickness".  In  order  to 
compensate  this  work  and  to  put  us  in  accord  with  reality,  we  will 
Introduce  a new  distribution  of  source  of  a certain  form,  which 


10 


will  bring  the  wing  back  to  a mean  nought  thickness.  The  variation 

of  vertical  velocities  w^q,  equal  with  -wjq^*  respectively 

at  the  extremity  of  the  wing  (where  the  indices  (1)  correspond  to 

the  two  marginal  parts  of  the  wing)  will  correspond  with  "the 

wing  compensating  for  slope"  of  symmetrical  thickness  which  will 

cancel  the  slope  and  the  effect  L of  the  slope  J — from  the 

X i 

extremity  of  the  wing  2).  In  this  way  the  term  which  produces 
a velocity  tending  towards  infinity  at  the  leading  edge  Is 
eliminated  according  to  a logrlthmlc  expression.  The  distribution 
of  sources  q>y  will  be  necessary  to  create,  at  the  same  time,  on 
the  lower  side  of  the  wing,  a distribution  of  pressure  without 
accentuated  peaks,  approximately  constant  with  the  exception  of 
the  regions  near  the  leading  edges. 

Then,  for  simplification,  we  will  chose  the  following 
functions  for  the  distribution  of  sources: 

5-20  = Arao  iq(1  -f-  *10  ^),  (20a) 

q'n  = Ki  *)(1  + *u  »))  (0  < T)  < l).  (20b) 

Thus,  we  got  those  two  "thick  wing"  components,  2).  and  3)*» 

In  order  to  be  reduced  to  one  with  the  slope  variable,  having 
therefore  the  mean  nought  slope. 


For  the  "compensating  wing",  the  function  of  the  axes  velocities 
Uic  will  be,  similar  with  (17),  the  following: 

nu = i *.  - ft** + *>  co,i>-f a d” + 


ii 


r 


+ T S„  “ ««  *>  C°8h_1  f(1  d>1  ~ X-  (21) 

In  the  course  of  the  calculations  we  are  driven  towards 
the  expression 

i OL.  = -^4-  f B* (l*  — x *)  (cosh-1 1 + cosh-1  + 

x,  ‘ 2kB*L  l B(l  — x)  B(l  + x)J 

+ 2 B*  x*  cosh"1  — )-  2 ^1  — ]/l  — B*  l*  j |^1  — B*  as*j  + 


+ 


h*h?\B*(i*— x*)  cosh-1  ^-^  + B*(I*+**)coHh-1^^  + 


3ttB*L 


B(l-  x) 


B(l+x ) 


+ (2  sin-1  Bl  B*  x*  + sin-1  Bl  — Bl  U — B*l*)  ^1  — B*  x*j  + 

+ s [«*(!*— x*)  fcosh-ll~-g--- cosh-1 1 + Bt~)  + 

\27rB1  L l B(l  — x)  B(l  + x)J 


+ 2sin-,Bt  Bx 


+ kn’h}\B»{l»—x»)  cosh-1  — — — B»(l*  + *»)cosh-1^^  + 


3ttB*L 


B(l  — x) 


B{l  + x ) 


+ 2 B*x*  cosh  - 1 + 2(l  — ^1  — B*ll)  Bx  \l  — — L.  (22) 


Through  the  superposltlonlng  of  the  three  wing  components, 
the  resulting  Imaginary  wing  Is  obtained,  equivalent  from  an 
aerodynamic  point  of  view  with  the  real  wing,  for  which  the 
axis  of  disturbance  velocity  will  be 

qix  = ruu  + ^ 


We  observe  that  the  velocity  Uj  j on  the  higher  side  of  the 
wing  Is  equal  and  of  opposed  sign  with  that  of  the  lower  side, 
corresponding  to  a thin  wing. 


12. 
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THE  SIMPLIFIED  CASE  OF  CONCENTRATED  SOURCES 


If  we  will  consider,  for  simplification,  a brief  variation 
of  vertical  velocities  in  the  abscissa  point  n « c due  to  the 
Influence  of  the  vortex  from  the  right  of  this  point,  we  will 
have  the  case  of  a brief  drop  of  Incidence,  corresponding  to 
an  "edge"  of  separation  of  vertical  velocities  Mo1  + < e)  x,  and 
(»!i’ + toll* c)  a,  fa J.  . The  brief  drop  In  the  distribution  of 
incidences  will  correspond  to  a variation  of  vertical  velocities 
expressed  analytically  by  the  concentrated  source  In  s»c  of 
Intensities  Qq  - <£>♦  , for  the  straight  part  of  the  wing, 

where  Qq  corresponds  to  natural  symmetries,  and  Qo  to  forced 
incidences. 


Following  in  continuation  the  path  presented  above,  as  In 
(1)  and  (3),  *e  will  obtain  for  the  three  wing  components  the 
following  expressions  of  the  axis  of  disturbance  velocities; 

1 « •>  1 2 


•x) 


x) 


_ ^io  + 2_ 


yirzgr-  - ~ [Qm  C08h"‘  ± 

±Qu*ooah ‘yj/^rzrsi]  (24) 

for  the  "lift  wing  ", 

i * - T««- + «•* »>  ^'£7)  + 
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+ f'  +a bVI+*)BX)  + l“ 

- ± [«„  oo»i>-  ]j  ± «„  * *« ) [ g-* !!::,] +i 

(25) 

for  the  "wing  of  symmetrical  thickness",  haring  the  slope 
equal  with  the  Incidence-^-  of  the  lift  wing,  which  was  deduced 
from  (17)  considering  the  sources  of  Intensities  Q20  and  Q2l, 

concentrated  In  jr  * c in  place  of  linlar  distributions,  (2),  and 

ai'u^aiu  (26) 

for  the  "compensating  wing"  of  slope,  where  ulc  has  the  expression 
given  by  (22). 

We  will  remarlcthat  given  as  the  sign  of  the  axes 
velocities  Is  different  on  the  two  sides  of  the  wing.  The 
expression  of  the  total  axis  vleocltles  will  be  that  given  by 
(23),  which  will  give  us  the  distribution  of  pressure  on  the 
deformed  wing. 


3.  THE  DETERMINATION  OF  THE  CONSTANTS 

For  the  determination  of  the  constants  q20  and  q21,  which 
appear  iu  the  expression  of  ^^(16),  will  begin  from  the  conditions 


Im  C(  Kl  — B*  x*  — f — iiA  da;  = ^!S , (27a) 

Jo  da;*V  d>)  ) d>) 

deduced  from  the  theory  of  conical  motion  (4),  the  Integration 
being  accomplished  on  a Semicircle  6 of  very  small  radius  around 


a certain  point  y ■ n oii  the  wing,  contained  In  the  Interval 
(y  ■ a,  y ■ 1)  (fig.  2). 


i -y +<» 


D,  h 


0 S A *, 

■rrmT 


obtain 

the  relations 

die;,! 

(28  a) 

drj  J’ 

’ail 

d* 

(28b) 

which  establishes  the  dependence  from  the  variation  of  sources 
and  the  distribution  of  vertical  velocities  on  the  thin  Imaginary 
wing. 


Starting  from  these  relations  and  keeping  in  mind  (l4f, 
we  place  the  conditions  at  the  limit  In  the  points  n « s and 
n ■ 1,  for  the  vertical  velocity  w' ■ wjq*i  ♦ woix2»  and  we 

obt*ln  «. 

I J,  r)  (2 B*  T)*) 


...  _ «...  _ Ml  ('  " — " V)  , 

" 1 ).  ,|2--S*V)  ’’ 

m _ J a - n)  (1  - B*  if)«  dr), 


In  the  case  of  natural  symmetries,  for  which  q^0  ■ 420  • and 
Wqj  * 0.  Here  we  find  the  first  relations  among  the  constants 

««T>  «i  < : X 


15. 


4 «&'"  [ V1  — B'  i*  — Vl  — B*  t*  — cosb"1  — + 

2 L Bl 

+ cosh"1  ~ — |l  — yj  V*  — -B*  «*  -r  4 Bl  — sin-1  Bi ) - 

_ 4 U-  ipM! + ^ i+MI_ain-. 

2 ' 1‘2  — Bl  Vi  + Bl  \i  — B» 

_sin-.L+£iO\  + ‘ L„-  . 

V2  + B#J  y2£il  \2  — Bl  Vi  + Bl 


:il+Bsf2\  1 /.  1 — BlVi 

\i  + B$)  + )'2  Bl  (8ID  Vi  — Bl  81 

_8in-,  IpM*  + sin*  4^1]  = 

V2  — B»  Vz  + Bt)} 


Vi  + Bl 


— ^ J^3  (sin"1  Bl  — sin1  6«  + 6i  Vl  — B*  l*—  Bt  Vl  — B*  «*)  + 

+ 2 (.61(1  — 6*  l*)*1* — B»  (1  — B*  s*)*11)  + — ((1  — 

5 Bl 

— ( i — b*  «*)  «*)J  = fc',y  — tc\y . (3ot 

For  the  Incidences  with  forced  symmetry,  making  q^O 

/ ff)  / 

q20  and  wjo  * 0»  from  (28a)  and  (28b)  we  obtain 

- s?  C ||-”>  d„  I» 

I J.  n 

■g— <■’—  (»> 

BH  J.  i)* 

from  which  we  deduce  the  following  relations: 

«y  *.  j|i  + -]  vr—  b*t* — 2 vi  — b'T* — 4 Bl  (8in_1  Bl  — 8in_l  Bt) — 

_ 4 6*  p (vr=w*  - f vr=w) — -y[(i— 

jja  s1)^*]  + cosh"1 4 — cosh"1  = a’ii*  — w®»'  * (32a) 


$'r,  (<•  - «■  + [*•  ‘'(3  '• r)  + ?(1-8  *)]*=*»  + 

4-  3 Bl  jBI^coBh-1  A_  C08  l^|  + 8in‘1  Bl~ 8in"IB8JJ  ■“ 

= tcJS'— «**»'.  (32b> 

Oh  the  other  the  mean  vertical  velocity  or  the  Incidence 

of  the  real  wing,  given  by  (3)»  equal  with  that  of  the  first 
wing  components  as  with  that  of  the  resulting  Imaginary  wings. 

Is  obtained  starting  from  relatlon(2).  In  this  way,  keeping 
In  mind  that  the  two  thick  wings  2).  and  3)*  are  compensated 
reciprocally,  the  relation  (2)  Is  written  simply 
4*  [ Mo  + w'n  *))  dT)  = *®io  + 4"  *» 

| Jo  " 

corresponding  only  to  thin  wing  component  1). 

Next,  the  relation  can  be  put  In  the  fora 

«C  » + **  + ( (»!.  -f  w«tq)  di)  = «>10  l + 4"tr®i  **»  <34> 

^ J«  z 

which,  after  accomplishing  the  Integral,  becomes 

.• , _ (A  [ B (61  - 5 »)  (1  - B*  «*r  - Bl  (1  - B * l*)*'*]  - 

— 2 [£2(1  — BH3)313  - B»  (1  — B*  a*)*'1]  — 

_3  (Bl  B*  l*  — B»  \ 1 — B3a3  — sin-1  Bl  + sin-1  Bs)  J * (35> 

In  the  case  of  natural  symmetries  of  the  Incidences  and 

tfoi'  — Woi  = ^S~|s  Cain”'  Bl  — sin'1  Bs  + Bl\l^B*l*  - Ba^i—BU3]  + 

+ 2[Bl(l  — B3l3)313—  Ba(l  — B3a3)313]  -f  ~[(i  _ B*l3)113  - 


— (1— 


for  the  case  of  Incidence  with  forced  symmetry 


Per  the  determination  of  the  constants  A^and  A 22  which 

appear  In  expression  (l6),  we  will  take  Into  consideration  the 

w' 

variation  of  vertical  vleocltles  ” ■ w^Q  + wQly,  on  the  first 

wing  component,  from  a point  on  the  wing  to  one  of  nought  vertical 
velocity  (for  example  on  Mach  cone),  as  Is  proceeded  In  the  theory 
of  conlc&l  motion.  As  In  (1)  and  (3)*  In  order  to  avoid  some 
difficult  calculations  to  determine  these  constants,  we  will 
consider,  through  approximation,  thsit  the  sources  distributed 
in  a line  In  the  interval  (s,/)  are  concentrated  In  y ■ b(  of 
Intensities  Q2q  and  Q2i»  8UCh  we  have 


(371 


Proceeding  in  this  way,  we  can  write  the  relations 

_,nr-“  tEDt*  - <■,  (w 


Im 


Juipk  * 

<*o#icuJ  (1  It) 


J.tU4  (I) 


dx* 
d •Wu 


h — B*  *»  — — ^-dx  = tril*, 
dx* 


(38b 


where  Un  represents  the  axis  of  disturbance  velocity  for  the 


simplified  case  of  concentrated  sources  In  y 
expression  (24). 


s^  given  by 


Through  the  accomplishment  of  the  integrals  which  appear 
above,  on  a circle  of  very  small  radius  around  the  origin  for 
(38a)  and  on  the  axis  of  the  abscissa  (^  - 0),  between  the  limits 
/and  1/B  for  (38b),  results 


18 


Aw  + 2l*An  + cosh-  ~]  + 

+ ?« 1 [2(c°8h'‘  7 “ F-?)  + f V1 ~ V 6081  f ]}  - 

-T"*"'  . (* 


A„-r^uP)^Zr^mk)~Klk)]  + A»EW- 

-4  (l  - T r F+T(4  2,8  ** y * + nw** **)  = 


(38b) 


In  which  we  noted 

it  *' 


r 


2 ftj*— «'*)  (1  — W») 


»»  #'*  1 

— — — n ( P,  *:)  + — 1 


b*«'» 


lfw, 


(40a) 


lll?,k|  - Kit)  + B(9.,  f)  + 

+ (K(k)  — E(k))F(9o,  fc')j,  (40b) 

(completed  ellptloalf 

and  K{k),  £(*)^ff  p (Pffc)  represent  theM.ntegral s of  the 

first,  second  and  third  instances  respectively,  having  the 
module  k and  parameter  p given  by  the  relations 


l =yi  — B*l\  P = B*s'*  — l = — B*(l  + 2*)*  — 1,  *'  = Bl,  9,  = 

9 

= 8in>i--  (41) 


Due  to  the  separation  of  flow  at  the  edges  of  the  wing 
and  resulting  vortexes  on  the  higher  side  of  the  wing,  fintte 
velocities  in  those  points  are  realized.  Imposing  these  conditions 
we  will  be  able  to  write  the  relation 

A»  + i*^t,  *=»  0.  (42) 


19. 


1 

Next  , it  is  observed  that  lnthe  expression  of  - U_  given 

am 

by  (22)  appears  the  constants  k1Q,  k 2Q,  kn  and  k21,  which 
follow  to  be  determined.  Taking  into  consideration  the  roll 
of  the  third  wing  component , which  will  have  the  mean  slope 

-*10*1  * woix2*  and  k««Plng  ln  «lnd  relations  (20a)  and  (20b), 
we  can  write,  similar  with  (33): 

— j (w»i  + M’oi  >i)  = — i ol% — v>u  1,  (43) 

from  where  we  deduce  the  relations 

(IT, o — »li‘)  l = - iSS'  ( V (1  +k\V  n)  (!—•«*  V)M  di,  (41* 

Jo 

pi  n R*  ri*)*11  . 

(«„  - »»)  I = *JS'  V (1  + 1>  r)(2_B,^  dT>’ 

in  the  case  of  incidence  with  natural  symmetry,  which,  ln  the 

course  of  accomplishing  the  Integrals,  becomes 

ir,0  — tcii>  = (1  — B*  l*)*1*  — JL(i_  Bii«)»/«_ 

224*  \ 3 12 

-±^r=TlU-  + _(1  _ _ 


--±-(l_(l-B«  !»)«•)] 


}■ 


(«2L 


- <!  = ^ [ Bl  V 1 - BH*  (3  + 221V)  - 11  si 


sin"1  Bl  — 


1^2  — Bl  • K2  + BlJ] 


+ 


+ — [B*i*y  1 — b*i*— 2(i— yi—sn*)]  + 

3 


+ 2n  — 4 ^sin'1 


1— Biy2 

V2 —Bl 


+ Bin*1 


l±BiK2\l 
V2  + Bl)  J 


(Tfe 


20. 


For  the  particular  case  of  forced  symmetry,  we  will  start 
similarly  from  the  relation  (42),  keeping  In  mind  and  wellknown 
from  (20a),  (20b),  (28a)  and  (28b),  and  obtain 

(wn—  u4V)  l*  = rj  (l  + KU)  (l  — ^ 


(«« -<’)  P-  — (1  + fcSti)  (l  — 

X)1  Jo 


which,  after  accomplishing  the  calculations,  terminates  In  the 


form 


u-iV — W0l 


ji  (( i —Bn*)**—  i ) + i.  o J-i  _ 


12 


-il' 


i_M*  + 


win  1 Bl 


ttiV—  W’oi  = -J§f  [y  Bi[a(  1 - BV)*'*  + 3 (if  1- BH*  + - 

_M((1_*.|.)*_l)j.  (47b) 

We  will  remark.  In  continuation,  that  equations  (30a), 

(30b),  ( 35 ) • (39a),  (39b)  and  (42)  constitute  the  system  for 

the  determination  of  the  constants  for  the  curved  thin  delta  wing 

(with  natural  symmetry).  Thus,  In  (39a),  in  which  we  replaced 
* *(r0  c«>) 

q for  q and  Wqi  « 0,  and  In  (42)  results 

a-  - ! “ [(* + f ) K r-l  r - - 7]  • + 

+9,T<  <«.) 

Iiitl 1 . ln> 

(48b) 

and  from  (35  and  (39b)  we  yeduce  the  relation 


1 


L 


Y—h'lk ) 


,+2t 


-2 (Bl  (l—BHi)^—Bn(l—Btiii)tlt)  - 


-■■lill/ 1 — BH* — if#  ^1  — if2#* — sin-1  Bl  -f-  sin  1 B»)  j | = _ » 1# ; (49) 


then,  together  with  the  relation 

4 Bq^l/l—BH*—  \iZlBh*—  cosh'1  — + couir*  — — 
l Bl  if# 

— ( 1 — y ) /i—  fiM  + ^ (bin'1  Bl  - siir  ‘if#)  — 

If.  _.l —Bl\2  , . ..l-fiff/a  . i—  B#f2 
— Ism  1 4- 8in  1 I — ain-i 1 — 


-sin 


2|”“*  \Z-Bl  ‘ ™ 1/2  + Bl  If2^—Bt 


. _.l  + fi#^21  , 1 f . l—BlV‘2  . 1 +B11/2 

— sin  1 — — ~ I 4-  — — I Kin  1 ’ sin  1 — L_ 


^2+ if* 


[■ 


Bin 


\2Bl  L”  |/2  — Bi 


Bin 


^2  + B/ 


/2 — if#  ^ K2+B*]} 


+ qT  J3^in-‘BI  — Bin”1  B»  + Bl  Kl  - B*l*—B»  Kl—B*«‘j  + 
+ 2[Bt(l—  BH*)*1* — B#(l  — B*«*)W]  + 


+ “ [(1  “ - (1  - BV)*»]J  = 0, 

5 in 


(* 


obtained  In  (30a)  and  (30b),  Is  formed  the  system  of  the  two 
equations  which  determine  the  constants  q^Q  AND  q*^ 


Proceeding  In  this  way,  we  will  find 

Iff 


qT  = - 


qh"  = 


iffjff— iff«iff 
iff 


Wi, 


iffjff— iff«iff 

22. 


«>10* 


(51i 

(filV 


where  we  made  the  notations 


Jfi»=  2B  J2 JV(1 — B*l»—  y)Kl—  &*•*—  COBh‘l  + coHh'1-^ 

1#  = 3(sm"»  Bl  — sin-1  Ba  + Bl  ^ 1 — B*I*  — Bs  Kl — B***)  + 

+ 2 [Bf(l— B*f  *)“'»— B*(l  — £*«*)•'*]  + 

+ JL [(l  _ EH*)'!*  — (1  — BV)6'2],  (Ki 

5 Bl 


^ !{((' +t)  -"--t]  «*' + 

+ i B2/*y(i  — ~J  ^2+jJ.  (52c)» 

:!i'A '! - (,  |-;'.'l ' '|':1' 1 

•RrFKtS1- 


. I A (£(«{  _ 5*)  (1  — BV)8'*—  Bl(l  — B'lT*) 
1 [ 5 


9UBais| 

— 2(Bi(l  — B*f*)8/J  ~ B*(l  ~ BV)*")  — 

__  3(  Bl  \'l  — B*l*  — Ba)fi  — B*** — sin-1  Bl  -f  sin"1  Bt)  j J.  (52dV 


Froan(35)*  (45a)  and  (45b),  considering  for  k10  and  k 
suitable  values, 

L _ L 

(53f 


* __  1 
U“  21 1 


results,  for  the  constants  k^Q  and  kjjj,  the  expressions: 

; .1  = 2*l!  I JL  [B(6l  — 5#)  (1  — B*»*)M  — Bl(l— B*ia)w]  — 

6BM5 

— 2[B((1  — Btlt),lt  — B»(  1 — BV)*'*]  — 

*3. 


— 3(Bl\'l  — BH*  — Bn  fl—  B***~  8in_lKl4-8in“,**)J  X 
x |1  j BH*  yV^BH*  — 2(1—  Kl— BH*)  + 

+ 2"-4(' 


. 1—  BlV2  . .,..,1  + _ 

B,n  V2—  Bl  + \2+Bl  1 


_ — (3  + 2BH*)  — 11  sin-1  HI]  ■ 

2 Bl 

. 1 _ . l+iHlTaU"1 

^2  — Bl  )/2  4-  Bl  I] 


(64a)» 


*21  = — t #(«*  — 5*)  ( 1 — ZW)‘/*  — Jii(  1 _ _ 

— 2[Bl(l~  — i?«(l  _ fiV)8'*]  — 

— 3(£/|/fr — Vi8!8  — — Zf8*8 — sin-1  Bl  4-  8in_,fi#)j  x 

xi~(l~  BH*)*1* — 1 ( 1 — BH*)*!*—  1 (-L  &m~' Bl  4-  Ki=£V)  4- 
l *.  12  8 \ZM 

+ W*[\  l1  “ (1  - y(i  - (1  - *W)  p • (54b 


Considering  now  the  torsloned  wing,  having  a forced 
symmetry  of  Incidences  In  reference  with  plane  OxjX^  (fig/  1), 
the  system  which  will  determine  the  constants  will  be  formed 
from  the  equations  (32a),  (32b),  (36)*  (39a),  (39b)  and  (42), 
In  which  q • q V)  and  w^  » o.  We  will  obtain,  successively, 

from  (39b)  and  from  (42): 

21 8 


A (J) 

20  — 


*£(*)( 

it  IS 

T* 


-fj72  + f(l  BH'YqW+  KJ£  (55a 

1 +2t 


24. 


(55b 


and  from  the  other  equations  result  the  values  of  the 
constants  qgo^AND  similar  with  (51*0  and  (51b),  In  which 

It  1 8 noted 


+ Ly  i —B* «*— 2 Kl  — BH*  + cobU~*  ~ — cosb-1^-  - 


— — B/|bI^1— B2J*—  j Vl— wj  + 3(nin~lBl  — sin'*  Us) j — 


Bl  Bs 

lBl — sin"1 


_ ± [(i  _ BH*) 3/2  — (l  — B2#2)"2] } • 
3 


)■ 


(56a 


IS”  = (1  — E2i2)‘«  + iljl  — 2y  + B2s2  [2—  y]]  Ki  — £*»*  + 

-1  3BJ  fain-1  Bl  — Bin"lBs  + Bl  |coBh-1  — cosh"1  . (5*1 

j£=  [sin ~'Bl  — Bin"'  Bs  + Bl  YT^Bil*—BsYl^B*»*]  + 

Bl 2 1 8 

+ - [ Bl  (1  — BH*)*1*  — B*(  1 — B2*2)2'2]  + 

4 

+ — [(1  — B2l2)4'2  — (1—  B2#2)4'2]  + 

5 Bl 

-f  Bl  jc08-,y  + cosh-1  - |l  + — j^l y-  ~ 


BVY 

3E(k) 


(66c) 


J" = {“  ~ m'r  + £ [‘  - 2 ~ + B’  '‘(a  - y)]  fl  -BV  + 

+ 3 Bl  Jain-1  Bl  — sin-*  B*  + Bl  ^cosh"1  — oo»h-1  y-jl|  + 


-f  cob"2  -j 2co8h"‘  — + ^2 yjj/l — Ji  — 


Y — K(k) 


Hi) 


E(k) 


K)7 ^ 


(66d) 


25. 


The  constants  kS«^  and  iff)  ^deduced  from  (47a)  and  (47b), 

20  21 

keeping  In  mind  (53)*  We  will  have 

*£  = — Bqtf'  j3[8in-‘  Bl  — sin-1  Ba  + Bl  ]fT^BH*  — Ba  /l— £*i»]  + 

+ 2[Bi(l — B2i*)*'*  — B»(l  — BV)*'-]  + 

+ [(l  — ~ ( 1_  * V)M ]}  * 

x{l(l- B‘J*r  +1(1-  B*i*)W  + 

+ I(^  + ^)-ir.  ,57.) 

***  “ MHqlU'  jafsiri-1  Bl — «in'1B*  +BI  ^1  — ~B\'* — Ba  /fI_B*«*]  + ’ 

+ 2 [ Bl(  1 — B*J*)*/»  — B«(l  — B*»*)M]  + 

+ gb1[(j  -fiV)W-(l  ~ *•«•)•»]  j X 
x | 7 [2(  1 ~ + 3 ^i~BV  + + 

+ l((l-B*i,)»/*_1)p.  (67b 

THE  CASE  OF  CONCENTRATED  SOURCES 

Stratlng  from  equations  (27»)  and  (27*>)»  we  will  find.  In 
In  the  case  of  concentrated  sources,  for  Q2Q  and  Q21  the  following 
expressions: 

Qio  - 1—  [(2  — BV)  + c« (Sh 


(1  — B^*)*'* 
1 


(1  — B*c*)  M 


[(2  — BW)  «’-<»’)  + c« 
[<’  - <’  + B*c*<<' -<>)], 


(581 


26. 


F 


and  from  (38a)  and  (38b),  In  which  we  will  Introduce  Uji  given 
by  (24),  we  will  obtain 

At0  4-  21' A, t +--%  (QM  + 2cQu)  = - (» 

n c 1 « 

U,o  + A „ i»)  (£  (*)  - * (*))  + Au  E (k) - 

(59b 

— ^ [«.« «*<>*  y + °Qu  (y~K  (*» = «w, 

s 

where  Y la  given  by  the  expression  (40a)  for  s ■ c. 

Using  relation  ( 33 ) • we  obtain 

j(«C  + ~°  »ii’j  + — y)  |«®io'  + y (j  4-  c ) irivj  = »,o  + y u>oll,  (9 

which,  together  with  equation  (58a),  (58b),  (59a),  (59b)  and 
(42),  forms  the  system  of  equations  from  which  we  deduce  the 
constants. 


we  will  obtain 


In  this  way, 

2 B ' c2 yi'  — c'[(  1 — B'c*)  n ( p,  k)  -f  K(lc)  — E (fc)] 


<■ 

«’tn 


2 Bh']fl' — c'[(\—B*c')  n (p,  k)-\-K  (k) — E(k)]—nc |l  — jj(l  — B»c*)v« 

(61a) 

= 1 — 


M>io’ 


w 


10 


KCj(l  — B*  c*)% 


2 B'c'j/l*— c*[(l—B'o')  n ( p.  At)  -f- K(k) — £’(*:)]— no  |l  — -ij  (i_ B'c')* 

(61b) 


27 


2 b*c*  yi *— c*  [( 1 — B' c»)  n ( P,  *)  + K (k)  — £(*)]— nc  (l— yj<l  -B*«')* 

(Hit) 


in  the  case  of  the  curved  delta  wing,  and 

C (1—2  B*c*)E(k)  + B*c*  K (k) ( 

Wn  (1  — — -^r  + (1—2  B,c*)E(le)  + B*ctK(k) 

(62  a) 

(1  — B* o')*'*  -f  |/l— £ [(1  — 2 B*c*)E(lc)  + B*c*  K(lc)] 

*»<  r * ____ 

*“  (l— (i  _£*„*)•/. + l/l  — ~ [(l-2B«c')i(k)  + BV£(i)] 

{ J ' (62b) 

= - — i*  U«  l - * ~-2  «'  - *&')  1 . (62c) 

4a  = — f < i - <|  ~-2Df,1°,>  « - *»>1  * (62d) 

7t  [ (1  — B'o')*'*  J 


<M8-« 


'01 


|l  — ^j(l— BV)»  + — ^[(1— 2Btc*)E(k)  + B*o*K(  1)] 

(62c) 


q&=—von 


B*c • 


|l— ^](1— BV)W  + ^1  — £ [(1  — 2B»C*)  U(Jfc)  + B*c»Jf(fc)] 

1 («2f) 


28. 


1 


for  the  torsioned  delta  wing. 


The  constants  k20,  k21  Is  deduced  Immediately  from  (45a), 
(45b)  (47a),  (47b),  taking  Into  consideration  ( 53 ) * ( 6lb) , (62a) 
(62b). 


4. 


THE  CALCULATION  OF  THE  DISTRIBUTION  OF  PRESSURE 
AND  AERODYNAMIC  CHARACTERISTICS 


We  have  shown  above  that  the  axis  of  disturbance  velocity 
on  the  real  wing  results  through  the  superposltlonlng  of  the 
three  Imaginary  wing  components,  obtaining  formula  (23)* 


For  the  calculation  of  the  distribution  of  pressure,  the 
total  axis  velocity  given  by  (3)  (fig*  3)  will  be  considered. 
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Moving  along  to  the  calculation  of  the  coefficient  of  lift 
of  theming  w^m 111  make  the  observation  that  the  wings  of 
symmetrical  thickness  do  not  give  lift  so  that  only  the  coefficient 
of  lift  of  the  "thin  lift  wing"  will  be  taken  Into  consideration* 


C. 
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Keeping  In  mlnd(l6),  we  will  obtain  the  following  expression 
of  the  coefficient  of  lift  of  the  wing* 
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+ T* 1 [(t  ~ 7 (> " * $))  (‘  - W-  -■*  t]|-  “ 


In  the  case  of  distributed  sources,  and 

= hf  + (^*®  + T ^‘jF1  “ 7*  ] ' (66) 

in  the  hypothesis  of  concentrated  sources  in  the  point  y ■ c 
on  the  wing* 
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In  order  to  define  the  parameter  j , which  enters  into  the 
expressions  above  and  which  determines  the  limits  of  the  die. 
trlbution  of  sources,  we  will  remark  first  that  the  position 
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of  the  maximum  distribution  of  pressures  coincides  with  the 
abscissa  y ■ c of  the  center  of  the  vortex  nucleus,  as  is 
ascertained  in  other  ways  from  experience.  But*  making  the 
calculation  on  the  basis  of  the  distribution  of  sources  (14)* 
it  is  ascertained  that  the  peak  of  the  depressions  on  the  higher 
side  of  the  wing  falls  approximately  at  half  the  distance  between 
the  center  of  gravity  of  the  intensities  of  the  sources  and 
the  abscissa  point  s*  on  which  we  keep  in  view  to  determine  them. 
In  this  way  we  can  deduce  the  relation  between  s and  c: 

e^8+~(l-s)  (-1  = 1,2  0,2}.  (67, 

Next*  for  the  determination  of  the  positions  of  the  vortar 
nuclei  centers*  we  will  observe  that*  at  curves  and  small 
torsions  of  the  wings,  these  can  be  considered  perfectly  plane 
and  parallel  with  the  undisturbed  flow  U^,  so  that  the  vortemes 
falling  off  at  the  leading  edges*  in  the  form  of  horns*  can  be 
considered  as  absent  from  the  higher  surface  of  the  wing.  This 
work  would  set  out,  however  small  the  intensity  of  the  vortexes 
would  be*  at  very  great  local  velocities*  incompatible  with  the 
reel  effects  of  the  separation  of  flow  at  the  edges.  In  order 
to  avoid  this  matter,  we  must  allow  that  with  curved  and  very 
small  torsion  of  the  wings,  accordingly  for  very  small  values 
of  the  parameter  (of10  ♦ O)  or  more  exactly,  when  Of  jo  ♦ 

J(X the  position  of  the  vortex  nucleus  will  be  c » Z» 
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For  greater  curree  (torsion)  on  the  wings  we  will  allow 
that  the  nature  of  vortexes,  which  start  on  the  leading  edges, 
evolve  proportionally  with  the  Incidence  corresponding  to 
respective  sections,  as  well  as  with  the  span  of  the  wing  (as 
with  the  wing  of  constant  Incidence  (1)),  accordingly  with  the 
square  of  the  span  of  the  wing,  and  winds  Itself  on  an  axis 
representing  the  line  of  the  centers  of  gravity,  which  correspond 
at  an  abscissa 

i.aj L.  («i 

l 4 

Ve  will  reaark  however  that  the  intensity  of  the  vortex 
nucleus  Is  something  smaller  than  the  total  vortex  generating 
Intensity,  the  rest  being  In  the  thin  layer  which  forns  the 
surface  of  the  horn  along  the  leading  edge. 

Fron  this,  directing  ourselves  according  to  experiments 
accomplished  on  the  delta  wing  with  constant  Incidence,  the 
minimum  position  of  the  vortex  center,  corresponding  to  curves 
(torsions)  or  large  parameter s 0( \ q can  considered 

o a 0,65  l -r  0,65  l.  (6D 

What  makes  this  minimum  position  very -stable  Is  the  fact 
that  at  great  incidences  Interior  vortex  nuclei  appear,  of  opposed 
direction  with  the  principal,  situated  between  this  and  the 
leading  edge. 
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Between  these  limits,  keeping  in  mind  experimental  results 
obtained  by  diverse  authors,  on  the  plane  delta  wing  with 
constant  incidence,  we  will  allow  for  the  position  of  the  vortex 
nucleus  the  following  approximate  formula  of  variation  with 
the  incidence: 


<» 


t “ l+M 

where  ^represents  the  Incidence  of  the  curved  wing  in  a suitable 
section  corresponding  to  the  center  of  gravity  of  the  aerodynamic 
effects,  namely 

3 (I! 

*~T 

For  the  torsloned  wing  (with  forced  symmetry),  due  to  the 
variation  of  incidences  in  section  x ^ in  function  x2,  we  will 
consider  like  the  incidence  calculated  in  the  formula  (70), 
lKP r*  point 

3 2 1 , 

X,  — ■ X,  = lx  i = —l- 

*4  3 2 


(72) 


We  will  remark  that  the  whole  rational  which  led  to  stability 
formula  (70)  can  be  applied  to  each  section  of  the  wing,  obtaining 
in  this  way  a curved  line  for  the  positions  of  the  vortex  nuclei. 
However,  deviations  from  a straight  line  are  ascertained  at  small 
incidences,  towards  the  peak  of  the  wing,  where  the  contribution 
of  lift  power  of  the  wing  and  the  distribution  of  pressures  is 
small  due  to  small  iurface  and  incidence,  and  in  the  back  of  the 
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wing,  whore  the  Incidences  grow,  the  vortexes  are  situated 

rectlllnlar.  In  this  wap  we  have  succeeded  In  planning  systematically 

(which; 

a vortex  sheet  in  the  fora  of  a horn,  on ’'an  axis  - being 

\ A'f  W 

considered  a straight  line  -Srthe  effects  on  the  wing  m the 
vortex  nuclei  are  found ^introduced  through  distributed  sources. 

This  work  was  necessary  because  only  in  this  manner  can  we 
apply  the  methods  of  conical  flow  of  the  higher  order  (2).  *n 
opposed  cases,  the  nuclei  being  situated  on  a curved  line,  a 
quaslconical  vortex  sheet  results,  and  therefore  motion  becomes, 
as  such,  more  complicated. 


Received  by  the  Editor  on  Dec.  10,  1968 
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THE  THIN  DELTA  WING,  WITH  VARIABLE  GEOMETRY,  OPTIMUM 
FOB  TWO  SUPERSONIC  CRUISING  SPEEDS*) 

by  ADRIANA  NASTASE 

Institute  of  Fluid  Mechanics 

Academy  of  the  Socialist  Republic  of  Romania 

In  the  present  article,  the  form  of  thin  component 

l with  variable  geometry^/ 

surfaces  of  the  heavy  lift  delta  wing^ls  determined, 
which  has  the  wing  span  2 1 given  when  it  is  completely 
folded  (fig.  1)  and  which  must  be  optimum  when  it  is 
placed  in  the  flow  characterized  through  the  MACH  number 
M,  and  when  it  is  completely  unfolded  (fig*  2)  it  has 
the  wing  span  2 L>  2 1 and  must  be  optimum  when  it  is 
placed  in  the  flow  characterized  through  the  Mach  number 
Mj<  M. 


*)  The  following  principals  presented  in  this  paper  were 
communicated  at  the  Mill tar ^Teohni cal  Academy  in  Brno.  Nov.  1968, 
under  the  title  "The  Thin  Delta  Wing,  with  Variable  Geometry, 
Optimum  for  two  Supersonic  Cruising  Speeds". 


The  designs  for  the  thin  delta  wing  with  variable 

geometry,  the  lift  powers  and  moments  of  dive  In  the 

folded  and  unfolded  position,  as  well  as  the  fora  of 

the  central  airfoil  are  given.  In  order  to  avoid  the 

falling  off  of  flow  on  the  subsonic  leading  edges  we 

will  cancel  the  axis  of  disturbance  velocity  u along 

the  subsonic  leading  edges  of  the  wing  In  the  folded 

and  unfolded  position.  Similarly,  the  wing  In  unfolded 

position  has  to  be  continued  at  the  traverse  of  the 

U^nfQldedJ 

edges  of  the  folded  wing  (becoming  edges  on  the* 
wing).  Through  the  application  of  variational  methods, 
the  problem  of  the  fora  of  the  surface  of  the  thin  delta 
wing  optimum  for  two  supersonic  cruising  speeds  Is  reduced 
to  the  resolving  of  an  algebraic  system  of  llnlar  equations. 

1.  GENERAL  CONSIDERATIONS.  THE  REDUCTION  OF  THE  PROBLEM  OF 
EXTREMUM  OF  THE  THIN  DELTA  WING  WITH  VARIABLE  GEOMETRY,  PROPOSFD 
BY  THE  STUDY  OF  TWO  VARIATIONAL  PROBLEMS  IN  CASCADE  REFER MG  TO 
THE  FOLDED  AND  UNFOLDED  WING 

From  Investigations  of  specialized  liturature  we  have  never 
found  another  theoretical  paper  referlng  to  the  determination 
of  the  form  of  the  surface  of  wings  with  variable  geometry, 
optimum  for  two  supersonic  cruising  speeds.  On  the  other  hand, 
wings  with  variable  geometry  present  a special  interest  In  modern 

3b. 
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aviation,  since  they  have  a series  of  advantages  over  the  wings 
•f  the  Invariable  fora,  advantages  among  which  we  mention  the 
most  Important: 


— facilitates  the  landing  reducing  the  length  of  the 
runway  and  the  landing  velocity; 

— permits  the  Increase  of  lift  capacity  of  the  wing  at 
great  velocities; 

— permits  the  control  of  aceleratlon  and  thermic  control 
exit  and  entrance  In  the  dense  atmosphere. 

We  will  study  In  this  article  the  form  of  the  surface  of 
thin  delta  wings,  with  variable  geometry,  optimum  for  two 
supersonic  cruising  speeds.  This  type  of  wing  can  be  found 
being  utilized  especially  In  aerospace  vessels  which  exit  and 
reenter  In  the  dense  atmosphere. 

In  the  present  work  we  begin  the  study  of  the  heavy  lift 
delta  wing  with  variable  geometry,  optimum  for  two  supersonic 
cruising  speeds,  eventually  provided  with  an  edge  of  central 
separation  for  which  are  given:  the  layout,  lift  powers,  the 
moments  of  dive  and  volumes  of  the  wing  In  the  folded  and  unfolded 
position,  the  form  of  the  central  airfoil,  as  well  as  a series 
of  geometric  conditions  enclosing  the  wing  along  its  leading  and 
back  edges,  continuity  of  height  of  the  unfolded  wing  along  the 
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leading  edges  of  the  folded  wing  ( becoming  edges  on  the 
folded  wing).  Similarly,  the  axes  of  disturbance  velocities 
must  be  nought  along  the  leading  edges  of  the  folded  and  unfolded 
wings  in  order  to  avoid  at  the  two  cruising  speeds  the  birth 
of  horn  shaped  vortexes,  which  have  the  tendency  to  form  along 
the  subsonic  leading  edges. 

Since  as  in  the  framework  of  llnlar  theory  the  effect  of 
lift  power  can  be  separated  from  the  effect  of  weight  in  the 
folded  position  as  well  as  in  the  unfolded  position,  we  will 
have  therefore  to  study  the  thin  component  and  that  of  symmetrical 
thickness  separately. 

The  study  of  the  thin  delta  wing  with  variable  geometry  is 
reduced  to  resolving  the  following  two  problems  of  extremum  in 
cascade:  first,  in  which  the  optimum  form  of  the  surface  of  the 
folded  wing  (of  wingspan  2 1)  (fig.  la,b)  is  determined  at  the 
cruising  speed  characterized  through  the  Mech  number  M,  and  the 
second,  in  which  the  optimum  form  of  the  added  surface  S (which 
is  presented,  in  our  case,  under  the  form  of  u>h.r|y  at  the  leading 
edge)  i 8 determined,  in  such  at  manner,  that  the  whole  unfolded 
delta  wing  ( of  wingspan  2 L)  (fig.  2a,b,c)  will  be  optimum  at 
the  cruising  speed  characterized  through  the  Bach  number  Mj. 

• . . - • • 
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The  first  variational  problem  consists  therefore  in  the 

determination  of  the  form  of  the  surface  of  a thin  delta  wing 

which  has  the  resistance  at  minimum  advance  at  the  flight  speed 

characterized  through  the  Mach  number  M and  for  which  are  given 

the  projection  in  the  plane  of  the  wing  (represented  through 

lsosclea  triangle  OAjA2 ) * the  lift  power  C]_,  the  moment  of  dive 

C_  and  the  central  airfoil.  In  addition,  we  will  consider  that 

the  axis  of  disturbance  velocity  u is  canceled  along  the  subsonic 

leading  edges  OAj  and  0A2  in  order  to  prevent,  at  the  considered 

cruising  speed,  the  birth  of  vortexes  in  the  form  of  horns.  In 

the  calculation  of  the  resistance  at  the  advance  we  will  also 

V^suctlonJ 

Include  the  forces  or^gmatlroTiB-  which  appear  on  the  subsonic  leading 
edges  of  folded  thin  wings. 

Also,  we  will  presuppose,  for  generality,  that  the  wing  is 
provided  with  a central  edge  OC. 

The  second  variational  problem  consists  in  determining 
the  form  of  the  added  surface  S (shown  in  fig.  2a,b,c),  such 
that  the  whole  unfolded  wing  will  have  the  minimum  resistance 
at  the  advance,  at  the  cruising  speed  characterized  through  the 
Mach  number  Mj<M  and  for  which  are  givens  the  design  of  the 
wing  (represented  by  lsosclees  triangle  OA^A') , lift  power  CL  add 
the  moment  of  dive  CM.  In  addition,  we  will  Impose  the  condition 
that,  along  line  OAj  and  0A2  (the  edges  of  the  folded  wing  which 


become  edges  on  the  unfolded  *lng),  the  height  of  the  surface 
of  the  wing  will  be  continuous,  and  the  axis  of  disturbance 
velocity  u will  be  canceled  along  the  subsonic  leading  edges 
OA^  and  OA^.  Similarly,  we  will  Include  the  advance  resistance 
and  the  forces  of  suction  that  appear  on  the  subsonic  leading 
edges  OA'I  and  OA^  of  the  unfolded  thin  wing. 


As  In  previous  works  (1)  - (6),  refering  to  the  determination 
of  optimum  forms  of  the  surfaces  of  thin  lift  power  systems , 
we  will  treat  the  two  above  mentioned  variational  problems  in 
the  following  two  hyposthesms  of  calculation! 

— the  first  hypothesis,  somewhat  classic,  in  which  the 
axis  of  disturbance  velocity  u is  infinite  on  the  subsonic 
leading  edges,  suction  appears  and  the  effect  of  the  falling 
off  at  the  leading  edges  is  omitted;  this  hypothesis  is  valid 
at  small  incidences; 

— the  seoond  hypothesis,  which,  at  the  chosen  flight 
speed,  we  avoid  the  falling  off  of  flow  on  the  subsonic  leading 
edges,  through  the  cancelation  of  the  axis  of  disturbance 
velocity  u along  that  edge.  In  this  case,  the  effect  of  conturlng 
of  the  subsonic  leading  edges  of  the  thin  wing  disappears  and 
once  again,  with  the  disappearance  of  this  effect,  the  forces 

of  suction  also  disappear  which  grew  ftreclsely  due  to  that  effect. 

As  has  been  shown  above  in  the  introduction  of  the  two 
variational  problems,  the  advance  resistances  and  the  suction 
are  also  Included  among  the  conditions  of  relation  of  the  problems, 
and  the  conditions  of  cancelation  of  the  axis  of  disturbance 
velocity  u along  the  subsonic  leading  edges. 

As  a result,  if  lnthe  final  systems,  which  give  the  solutions 

yx 


of  the  two  problems  of  extremum  mentioned  above,  we  exclude 
the  conditions  of  cancelation  of  the  axes  of  disturbance 
velocity  on  the  subsonic  leading  edges,  we  place  ourselves  In 
the  position  of  the  first  hypothesis  of  calculation*  If 
however  we  exclude  from  the  final  systems  of  terms  resulting 
from  suction,  we  will  place  ourselves  in  the  second  hypothesis 
of  calculation* 

2.  THE  EXPRESSIONS  OF  THE  AXIS  OF  DISTURBANCE  VELOCITY 
u ON  THE  FOLDED  AND  UNFOLDED  WING 

We  will  consider  the  wing  referred  to  a trlorthogonal 

system  of  axes  Ox^x2X3,  which  has  axis  OXj  parallel  with  the 

direction  of  flow  from  infinity  Uo  and  the  wing  is  considered 

plotted  in  the  plane  Ox-Xg*  We  will  presuppose  that  the  vertical 
(velocltleBJ 

*~'yof  disturbance  w and  w on  the  wing  are  expressed  in 
the  form  of  superimposing  homogeneous  pollnomials  wB-1  and 
w^j  in  Xj  and  x2»  that  is: 

N N m-1  f X,\  „ 

to  = £ to._l  = £ ar*  £ Wm_k_x.t  |y|*  It/  = •—  (H 

n.l  M«l  k — 0 V 1 ' 

on  OA.A.  and 

1 ^ N N m-I 

to  = £ U)._,  = £ asr*  £ |y  1*  (1 

»-l  »-l  1-0 

on  OA^A'  and  OAj^Aj  • In  this  formula,  and  wB_k-i.k 

ari  the  constants  and  N.l  lsthe  rank  of  the  pollnomlal  of 
highest  order  taken  into  cgnslderation  in  the  expressions  w and  7* 
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If  we  Integrate  expressions  (1)  and  (2)  In  reference  with 
we  obtain  heights  Z(xltx2)  and  Z(xj,x2)  on  OAjA2  and 
Oh 2^2  respectively  and  OA^A^  : 

Z(xt,  xt)  ^ — ■ j wdxl  + f(xt)  = j 

(s“f)  (* 


uwe 


Z(*„  xj  = — -Iwdx,  + f(Xt)  = 

V m J 


In  the  above  formulae*  f(x2)  and  ?(x2)  are  some  arbitrary 
functions  of  x2  which  we  can  eventually  use  In  the  placing 
of  supplementary  geometric  conditions  on  the  surface  of  the  wing. 

Considering  first  the  folded  wing  (closely)  (fig.  1)  having 
OA^  and  0A2  as  leading  edges  and  OC  as  central  edges,  the 
expression  of  the  axis  of  disturbance  velocity  u on  the  folded 
wing  will  be  (7) 
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for  the  folded  wing  with  both  supersonic  leading  edges. 


We  will  now  consider  the  unfolded  wing  (fig.  2).  Since 
the  geometry  of  the  central  portion  of  the  wing  Is  Invariable, 
It  results  that  w»  will  be  Invariable  on  portion  OA^A^,. 

The  former  leading  edges  OA^  and  0A2  become  edges  on  the  new 
unfolded  wing,  and  the  new  leading  edges  are  now  OAj  and  OA^. 
With  these  observations,  the  velocity  of  disturbance  u on  the 
unfolded  wing  will  be  of  the  fora 
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If  the  edges  OA'  and  OA^  are  subsonic  or 
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If  the  lending  edges  OA'  and  OA2  are  supersonic. 


Observations: 

a)  In  formulae  (6)  and  (8)  the  first  term  has  enlarged 

Itself  with  ^[T-B2y2  and  respectively  In  order  for  the 

radical  to  appear  In  the  denominator  and  so  that  we  can  bm.able 
to  better  systematize  the  calculation  of  the  aerdynamlc 
characteristics  of  the  wing. 

b)  If  the  edges  of  the  folded  wing*  as  well  as  those  of 
the  unfolded  wing  are  supersonic  for  the  Mach  number  of  flight 

then  In  terms  which  contain  the  coefficient  G^q  of  formula 
(8)  and  which  represent  the  contribution  of  the  edges  of  the 
folded  wing  In  the  expression  u "argch"  Is  replaced  with  "argcos". 
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Coefficients  An.2q,  Cn>2q  etc.  of  the  axes  of  disturbance 
velocities  u are  related  to  coefficients  of  the  vertical 
velocities  of  disturbance  through  the  linlar  relations  deduced 
from  the  compatibility  conditions  of  P.  Germain  (8).  Coefficients 
of  the  axes  of  disturbance  velocities  u on  the  folded  wing 
will  depend  only  on  coefficients  °f  wt  since  they  are 

functions  only  of  the  form  of  the  surface  of  the  folded  wing. 
Therefore  we  can  write 


(9a) 

(9b 

(*i 


Instead,  coefficients  of  the  axis  velocity  of  the  unfolded 
wing  will  depend  on  coefficients  of  w,  as  well  as,  w.  Therefore 
we  will  have 


n - 1 


2q  * Jl  (C ' 
*-0 

1 

(10a) 

C'n.2,  « 

o 

1-  CsT. 

(10b) 

- 1 

j-0 

g,T 

(10c) 

f-0 

1 h,T 

(10(1) 
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Observations 


a)  In  this  last  expression,  « 0 since  the  constants 

r 

Hnq  obtained  calculating  a semiresidua  around  point  y ■ L,  thus 
only  depend  on  the  drop  of  vertical  velocity  of  disturbance 
from  this  point,  which  is  -9.  For  homogenlty  we  have  written 
therefore  the  expressions  of  the  coefficients  under  the 
form  (10d). 

b)  In  formulae  (9a),  (9b)  and  (9c),  coefficients 
depend  only  on  the  form  in  the  plane  of  the  folded  wing  (thus 
of  1)  and  on  the  chdsen  flight  velocity  characterized  by  the 
Mach  number  M at  which  the  wing  is  folded  (B*')(m^-1  ) • 

c)  Coefficients  a^"1,,  a£\  etc.  depend  on  the  form  in  the 

plane  of  the  folded  wing,  as  well  as  the  form  in  the  plane  of 
the  unfolded  wing,  thus  by  1>  and  L,  as  well  as  on  the  Mach 
number  of  flight  M^t  which  the  wing  is  in  the  unfolded  position 
(B-tf^i.  Therefore,  for  same  designs  given  and  for  some 
flight  velocities  given  for  the  folded  wing,  as  well  as  for  the 
unfolded  wing,  coefficients  5*,"’/  etc.  are  constants 

in  the  study  of  the  two  problems  of  extremum  proposed.  As  a 

result,  from  the  first  variational  problem,  referlng  to  the 

folded  wing,  we  will  determine  the  optimum  values  of  the  coefficients 

w . . . of  w,  and  from  the  second  variational  problem,  referlng 
n-j-l . J 

to  the  unfolded  wing,  we  will  determine  the  optimum  values  of 
the  coefficients  wn-J-l.J  of  the  vertical  velocity  of  disturbance 
w.  In  the  second  variational  problem,  coefficients 

n 


L. 


of  the  vertical  velocity  of  disturbance  w,  determined  by  the 
first  variational  problem,  will  be  reckoned  constant. 

We  will  press  on  now  to  solve  successively  the  two  above 
mentioned  variational  problems.  We  will  consider  the  geometric 
conditions  of  relation  and  aerodynamic  characteristics  only  on 
half  the  wing,  due  to  the  symmetry  of  the  wings. 

3.  THE  STUDY  OF  THE  FIRST  VARIATIONAL  PROBLEMS  REFERING 
TO  THE  FOLDED  WING 


Before  moving  along  to  the  solution  of  the  first  variational 
problems,  we  will  explain  the  geometric  conditions  of  relation 
and  the  aerodynamic  characteristics  that  enter  In  the  problem 
of  extremum  of  the  folded  delta  wing.  In  this  direction  we  will 
use  a series  of  results  obtained  before  for  the  thin  delta  wing 
of  minimum  resistance  provided  with  a central  edge. 

3.1  EXPLANATION  OF  THE  GEOMETRIC  CONDITIONS  OF  RELATION 
OF  THE  FOLDED  WING 

The  geometric  conditions  of  relation  for  the  folded  wing  are 

a)  the  condition  that  the  central  airfoil  of  the  wing 
be  Imposed} 

b)  the  condition  that  the  axis  of  disturbance  velocity  be 

S'O. 
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canceled  along  the  subsonic  leading  edges  of  the  thin  wing. 

a)  From  the  condition  that  the  central  airfoil  be  Imposed 
(Z( x^ »0)«f(  ) ) through  the  development  lnthe  Mac  Laurln  series 

of  '■fMxj)  and  Identification  with  Zfx^O)  we  obtain  the  relations 

/(0)  =9(0),  (1U, 

= (llbj 

t\  - 0 - <r  (0)  = 0,  t ^ 1 , a, . . . n.  (1  lc) 

The  first  relation  (11a)  gives  us  the  value  of  f(0),  the 
second  (lib)  Indicates  that  the  central  airfoil  has  to  be  given 
In  the  form  of  a polynomial  of  rank  N,  and  the  relations  (11c) 
enter  as  relations  of  connection  In  the  studied  problem  of 
extremum. 

From  the  condition  of  canceling  the  axis  of  disturbance 
velocity  u along  the  aubsonlc  leading  edge  0A^( y*£,u*0)  we 
obtain,  similar  with  (6),  the  relations 


J-T  = £ = o, 


>«0 


*<;) 


a, 


x P41./,  1. -.•••**’■ 


(12a) 


Observations. 

To  assure  that  the  symmetry  of  the  surface  of  the  folded 
wing  In  reference  with  plane  OxjX^  and  as  a result  continuity 
of  the  height  of  the  surface  at  the  traverse  of  the  edge  from 

SV, 
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the  origin,  It  Is  encessary  here  that  the  function  f(x2), 
which  enters  into  expression  (3)  of  the  height  of  the  surface, 
be  parried  In  reference  with  x2,  In  other  words  to  have 

/( J'a)  - /(—  •'«)•  (13) 

3.2.  EXPA I NATION  OF  THE  AERODYNAMIC  CHARACTERISTICS  OF 
THE  FOLDED  WING 


The  aerodynamic  characteristics  which  enter  In  the  problem 
of  extremum  of  the  folded  wing  ares 

a)  the  lift  power  Cj; 

b)  the  moment  of  dive  Cn; 

c)  the  wave  resistance  Cd{ 

d)  the  coefficient  of  suction <5 q. 

a)  The  condition  of  the  lift  power  of  the  wing  to  be  given 
Ci  * const  becomes 

<\  E S'V.» ,'0,ls,•  (14> 

n I/O 

In  this  formula  are  constants  In  the  studied  problem 
of  extremum  and  which  are  of  the  form 


a;, 


*<J) 

£ 1) 


irmo>  \ D 1 ; <> 

* « i i ] 

Q 1 “7  ' 1 


(15) 


If  the  delta  wing  has  subsonic  leading  edges  and,  respectively. 
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if  the  delta  wing  has  supersonic  leading  edges. 

b)  the  condition  of  the  moment  of  dire  of  the  wing  to 
be  given  becomes 

N n-1  W )i  -I  1 , . 

r-  = £ £ r- «•.-»-> ■»-  £ £ )ITVV',,'~’  *•»  ' 

“o  n .1  ) »0  w f - 

In  this  formula  have  meaning  from  (14). 

c)  The  expression  of  the  coefficient  of  wave  resistance 


becomes 


.V  S n-lw-l 

r-=  £ £ £ £ iu  , >’ 

n J m ■ 1 i 0 k 0 


where£lninkj  are  the  following  constants* 


u„..,  = ^.+«.(o,o  + 

n / (Hi  + H)  I ,to 


Bl  h 


for  the  folded  delta  wing  with  subsonic  leading  edges  and. 


respectively. 
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for  the  folded  delta  wing  with  suoer sonic  edges. 

In  formulae  (15),  (16),  (19)  and  (20)  the  Integrals 
^7  ic  and  3k  ^axe  the  following  tralues 
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d)  In  conclusion,  similar  with  (1),  (9),  the  expression  of  the 
coefficient  of  suction  will  be  of  the  form 

N N tn  - 1 w 1 _ .1 

' - S S E E 
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3.3  DETERMINATION  OF  THE  OPTIMUM  VALUE  OFTHE  VEBTICAL 
VELOCITY  OF  DISTURBANCE  w FOR  THE  FOLDED  WINO 


We  have  accordingly  found 

Cd  m minim 

with  the  condition  of  connection 

c,  = c,„  = E,  — o,  F'v  = o, 


(25) 

t = 1 ,2  . . . A’ ; (2fi) 


similar  with  (1)  - (6),  the  calculation  of  this  extremum 

with  connections  is  reduced  to  the  calculation  of  the  extremum 


without  connections  of  the  expression 

//=('’„  4 r.)  -i  v"  (\  i-  x*'Cm  f-  £ (>,/:,  )-  (I, /•'!'’).  (27) 

t I 

Canceling,  as  In  (1)  - (6),  coefficients  of  the  variations 
<hw0«-»  we  obtain  the  following  equations: 

A’  n - 1 

y,  y (On. A K(+  1,0,1  I Oo+o*  I.H.jo  4 d»,n +o+l,n./  1“ 

n — I 1 0 

4"  do  KO  + l.nj.o)  M’»  1-1,1  4 l1'  ’A#  + 0 + 1,0  4"  I Ill'll  l.i  4- 
4-  n So  4-  iAfl  i o 1 1 Bo  i o f i ,o 

0,  IT  = 0,1... (A- 1),  I < fl  4 n I 1 < A.  (28) 


These  aquations,  together  with  the  conditions  of  connection 
(26),  form  a llniar  algebraic  system  which  determines  in  a 
certain  way  the  values  of  coefficients  Wg6  of  w,  as  well  as 
the  values  of  theTmultiplicators  and  u^,. 

The  first  variational  problem  referlng  to  the  folded  wing 
can  therefore  be  considered  solved* 


TW  3 

4.  THE  STUDY  OP  THE  VARIATIONAL  PROBLEM REFERINO 

TO  THE  UNFOLDED  DELTA  WING 


As  in  the  case  of  the  first  variational  problem,  we  will 
explain  first  the  geometric  conditions  of  connection  and  the 
aerodynamic  characteristics  that  enter  in  the  formulation  of 
the  problem  of  extremum  for  the  unfolded  wing.  After  that  we 
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will  consider  thesutfsmd  portion  of  the  wing,  that  is,  the 
folded  delta  wing  k\Okz  given  by  the  first  variational  problem, 
and  we  will  determine  the  form  of  the  added  surface  composed  of 
two  by  the  leading  edge  (formed  from  triangles  AjOA^  and 

A2OA2),  in  such  a way  that  the  whole  unfolded  wing  will  be  at 
minimum  resistance  at  the  second  cruising  speed  characterized 
by  the  Mach  number  . 


4.1  THE  EXPLANATION  OF  THE  GEOMETRIC  CONDITIONS  OF  CONNECTION 


OF  THE  FOLDED  WING 


In  the  case  of  the  unfolded  delta  wing,  we  will  impose 
the  following  conditions  of  connection  of  natural  geometry} 

a)  the  condition  of  continuity  of  height  of  the  surface 
of  the  wing  at  the  traverse  of  the  edges  OAj  and  0A2  (which 
are  the  leading  edges  of  the  folded  wing) ; 

b)  the  conditions  of  cancelation  of  the  axis  of  disturbance 
velocity  u along  the  subsonic  leading  edges  OAj  and  OA^  of  the 
unfolded  wing. 

Because  of  the  symmetry  of  the  wing,  the  above  mentioned 
conditions  will  be  put  only  on  edge  OAj  and  on  the  edge  OA^. 

a)  From  the  condition  of  continuity  of  height  of  the  surface 
\ unfolded) 

of  tKeWilhg  along  edge  OAj  (y«l,Z«2)  is  obtained  the  relation 

SL 


— w, 


>.,)  (a* 


which  is  similar  with  that  established  earlier  (10)  (in  which 
we  had  to  substitute  s »£)• 


If  we  presuppose,  In  a similar  way,  that  f(x^)  and  f{x^) 
are  developed  In  the  Mac  Laurln  series,  we  obtain  the  relations: 


(aoa 

m m v(*»,  (3oi. 

£ , |fc—  --  - 77  — /'*'(«)],  < - 1,2,.  • A'  (30* 

I / ® A o l — /*■  « • 


The  relations  (30a)  show  that  the  development  In  series  of 
the  arbitrary  functions  T(x^l)  and  f(x^)  have  to  have  the  same 
remainder  as  order  N. 

From  (30b)  results  that  these  functions  must  have  the  free 
term  the  same  and  equal  with  'sp(O),  as  well  as  the  results  from 
(30b)  and  (lla),  and  the  relations  (30c)  can  be  interpreted  as 
relations  of  connection  between  arbitrary  functions  t(x^)  and 
?(x2)  and  therefore  can  be  eliminated  among  the  conditions  of 
connection  oil.-,  the  studied  problem  of  extremum. 


b)  The  condition  of  canceling  the  axis  of  disturbance 

velocity  u along  the  subsonic  leading  edge  OA^(y«L,  u«0)  lead! 

to  the  following  relations  of  connection: 

£1. 

L _ 


«C) 

£ o. 


w = 1,2, 


If  we  mark  with  and  the  expressions 


«(i) 

b'  £ 


A ( !.  ) . 

e;,  - £ (3-b 

y U 

then  the  relations  (24)  can  be  written  under  the  fora 

• i i 

Fi"  = £ Kwv-^.M-fc), I - 1,2,  ...,Ar.  (■») 
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4.2.  THE  EXPLAINATION  OP  THE  EXPRESSIONS  OFTHEVCHAEACTERISTICS 
FOR  THE  UNFOLDED  WING 


As  in  the  case  of  the  folded  wing,  we  have  for  explainatlon 
the  following  aerodynamic  characterlsltcs: 

a)  the  lift  power  CL? 

b)  the  moment  of  dive  CM; 

n 

c)  the  advance  resistance  C ; 

D 

d)  the  suction  <S  . 


a)  The  condition  of  the  lift  power  of  the  wing  C.  to  be  given 
(C^  * CLq  * const)  becomes,  in  the  case  of  the  unfolded  wing, 
of  the  form 


C,_  =- 


^ ' w.r,  ri.r,  (\y  — 


= £ j)  (AJ, 4- A^fr.-j-u)  =»Ct,  = vonM- 


f (have  the  values  given J 

In  this  formula,  the  constants  \ Vin  the  . , 

r/iarej  *0  *tf*~*.t 

and  the  constants  ^^^obtalned  from/|'nj  replacing  the  coefficients 

ck  , c.'  etc.  their  corresponded ejbarrd.  Itfari. 
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b)  In  a similar  way,  the  condition  of  the  moment  of  dive 
to  be  given,  CM  * ■ const,  besomes,  In  the  case  of  the  unfolded 

wing,  of  the  form 


C„  = — — 

l (30  * J 


( « 

J<xm; 


rfcl.r,  (\>i  --- 


= 5j  Yi  ( r»/?r» = C„,  = const. 


_ / szj  ✓ 

In  this  expression,  the  constants  and  I are  connected 
by  the  constants  A and  /\n^  from  formula  (27)  through  the 


relations: 


I’'  w ^ 1 a ' r'  _ w "I*  1 . > 

11+2  v + 2 


c)  The  expression  of  the  coefficient  of  wave  resistance 
becomes.  In  the  case  of  the  unfolded  wing,  of  the  form 

r»=  y..t  I [ *c«.r,  da1,  fly  4 ( tCicr,  d.r.  d»/l  — 

U x1j  I J0A,C  ■ J 


K N m l n - 1 
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The  constants^/  and  l L . . are  given  In  , and  the 

— nmkj_  j n®*tJ  _ rxO 

constantsiinmkj  andJ/^j^^  are  deduced  froml/^^,  froa^/^^ 

respectively,  replacing  the  coefficients^,'01!  / ^ etc.  with 

i 2^,*j 


their  barred  correspondences 


d)  The  coefficient  of  suction  <T  is,  similar  with  (9),  (1), 

«•  - - ut  i i . ' ..  55  -»"■ — • w 

and  can  be  written  under  the  form  (1) 
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and  the  constants  a:.!^,,  a «r  are  obta?  '*1  replacing 

in  A,.,,,  , in  turn,  one  of  the  coefficients  ±,X>  or 

respectively,  both  coefficients  their  correspondences. 


4.3  SOLUTION  OF  THE  TWO  VARIATIONAL  PROBLEMS  FOR  THE  UNFOLDED 
WING 

We  have  now  explained  all  the  terms  which  enter  into  the 
problem  of  extremum  of  the  unfolded  wing.  Coefficients 
etc.  are  functions  of  the  design  of  the  folded  and  unfolded  wings, 
that  is  1 and  L,  as  well  as  Mach  number  Mj  of  the  second  cruising 
speed.  Thus,  these  terms  are  constants  in  the  studied  problem 
of  extremum. 


Similarly,  coefficients  w , , , of  w , which  have  been 

Tl-  J*  1 a J 

determined  from  the  first  variational  problem  refering  to  the 
folded  wing,  are  constants  in  the  second  problem  of  extremum. 


Therefore,  It  results  that  we  have  to  determine  the 
coefficients  of  " ln  8uch  a "ay  that,  for  the  unfolded 

wing,  the  expression 

V u | f s'  = minim, 

at  the  cruising  speed  defined  by  Mach  number  M1#  with  the 
conditions  of  connection 

CL  = const,  Cu  = const , 7<’l"  = 0,  t = 1,  2,  . . .,  A\  (42) 

This  problem  of  extremum  with  connections  Is  reduced  tto  the 
calculation  of  the  extremum  without  connections  of  expression 

li  = (Cu  \ o')  + Vt'j.  1-  >.u,<’w  + EM’"**  (-13) 


We  have  explained  all  the  terms  which  enter  Into  the 
hamiltonlan  expression.  Calculating  Its  first  variation  and 
canceling  It,  we  obtain  the  extrmum  of  H.  If  we  then  Keep  In 
mind  that  the  coefficients  w*a  are  Independent,  we  will  obtain, 
through  the  cancdlatlon  of  the  variations  S^tfcr  > equations  of 
the  form 


v;  % [(ni'ina+i.-M  + ««J .1 4 Ai>i + 

• 1 j-  0 

+ + “•+«.+ + Al'.l,,  + A04OhwJ  j t j + 

+ X‘»A„.+1-  + v*Tv0)l.0  + X0,a+1(-),!oi,.o  = 0,  (M) 

(0,  a)  = 0,  1 1), 

1 < 0 + o + 3 < JV- 


These  equations,  together  with  the  relations  of  connection  (35). 
make  up  a llnlar  system  of  algebraic  equations,  which  determine 
In  a certain  way  the  optimum  values  of  the  coefficients 
of  w,  as  well  as  the  values  of  the  multlpllcators  >1*  X+ 


d 


The  problem  of  determining  the  form  of  the  surface  of  the 
unfolded  delta  wlng;  optimum  at  the  cruising  speed,  characterized 
by  Mach  mumber  M^,  can  be  thus  considered  solved. 


5.4  Calculation  Appendix 

We  will  mark  with  Ijj  the  following  Indefinite  Integrals 

1 _ c ?/*  dy  . = _ ?*—  y wiOs  — ?/T)  + “-1  = 

' 1 1/  / ra 


, _ f _JL“| _ — v li\  (L2  — y 2)  + — 

‘ J V'«,  — ?/2)  k 

_ y (—  D*(fe!)  (2j)  [(^.ynr<.Riny  f 

h (j  !)3(A — J) ! /#,  L 

yfifcx*— »»)  f [(j— v) n*  ^ ' {T/  _ v),  vl 
+ _ .^l2v(2j_2v4  l)!v  ' J 


and  respectively  with  and  2;**  («,  p)  the  following 

definite  Integrals  of  type  I^s 

fP  yk  Ay  (46) 

71' = j,  y'  /?;-  ( l1  — j/2) 

(p  j/*cty  (47) 

These  definite  Integrals  can  be  obtained  taking  the 

Indefinite  Integral  given  In  (45)  between  the  limits  ^ 

and  making,  for  the  the  Integral^  the  substitution 


With  these  notations,  the  expressions/)  from  formula  (34) 

nj 

of  the  coefficient  of  lift  power  C will  be  of  the  form 

L 

a 1 * ( 2 ) *’(  2!)  r''"> 

- „.iir+-.i{ .5  M 4 "• L>  + 

+ SSHffi  s t a+tuir  'a;(„, 

45  «>ll  l»U  *1  T ' J 


If  fche  unfolded  wing  has  both  supersonic  leading  edges  and. 


respectively. 


....  « M+  CJ±L—  g:”  lo,  1 ) + 

“ v.Li.A  ill  A " r «,)  h -1  + 1 l *J 

4.  £ t 3;'-7M  t (-ini'  ‘ark  )}  (,,J1 

If  the  unfolded  wing  has  both  supersonic  leading  edges  and  the 
edges  (provided  from  the  edges  of  the  folded  wing)  are  subsonic. 


If  the  edges  OA  and  OA  become  supersonic,  then  In  the 

1 ^ _ _ 


expression  (46)\jl-B^i  is  substituted  with  in  terms 
which  contain  coefficients  g^n^  and  which  represent  only  the 
contribution  of  these  edges  In  the  expression  of  /1aj« 


In  a similar  way,  the  expressions^!  . and  , , from 

nmkj  nmkj 

formula  (3?)  of  the  coefficient  of  wave  resistance  of  the  unfolded 
wing  will  be,  respectively,  of  the  form 


Q „„  a. 

' VlL(m 


— J r a£'f  g[ , u (o,  i) 

»)  1 <i-o 


H-*1) 


• * ' /'<»>  1 f L*  ' 
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*-H II, k I 
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gy> 
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(Bll 


” ‘3:(l,  L) 


for  the  thin  unfolded  delta  wing  with  both  subsonic  leading 
edges. 


If  the  unfolded  wing  has  supersonic  leading  edges  and  the 

**  *0  **<1) 
and  ^ ^ i. « 

nmkj 

respectively,  of  the  form 


subsonic  edges,  we  will  have  the  expressions Q.  ”7  an<* nmkJ  » 


1 »1  nrecos  , W'  f)  jj(52 


cw\Ji 


!>*•'«  _ 
**•••»  I 


I ») 

*0f!) 


*C>  / 1 , 

S ('■ «,) 


' 2 > ,'im  / ( , [ i \ 

+ .5  A niP\ 


c '■ 


+ E- — f K i — «?  /=•  Vf  i 1 t (—  n*  'it* 4 '3r(i,  ') 

+ 2(-ir*f,  +(— i)*"i^""«w'iipl  ' 

8 /K( s) 

1C.,  = — t——  E ' 

V*mHm  ; »)  | 

R(  2 ) „>(m  / I f 1 ' N 


6V. 


, V \ k'  r="irf  If  V 1 1 4 (—  I )* " 1 11 " ' 1 X"  (0. 1 ) 

, i)  - (A  f </  4 1)1  < « 

4-  fl  4-  ( — 1 )*'  H * ' ' argi-li  [/  /{]  p 1 

4 2 ( — 1 )’  [ 1 4- (—I)1  "W*  ’ 1 :i>4M'li  |/  (l  u!‘\,)  \ f 


, 1 1 /)2  Ni 

4 2 ( — 1 )’  [ 1 4-  ( — 1 )*  "1  ’ 1 arjrrh  / , f.  , I 


+ V !SlL - ^ f Kaf  p -iicii  i ( — i >*  * ) v * * (o.  />  • 

f"«2(*-4  7 4-  1)\  « » 


4-  [1  4-  ( — IT’1]//  * 1 arccos 


1(1*  ' i — l"'  ')  iirocoR 


/()  | /*,  /-)(!  — «,M 
2 /{,  (/,  — /) 


+ v — isL-  *£  [1  +( — <^***fl»  i ) 

„Zj„2(Jt  4-  9 4-  D\  <-o  V *i/ 


' j5.u-.-o M + 

1 >• ' 1 1 1 — ( — i.)1-’1] areco. f ~I3B^1'++I)'i|l>  j}' 


If  the  edges  and  the  leading  edgeB  of  the  unfolded  wing 
are  supersonic,  we  will  have,  respectively,  the  expressions: 

ii,::.,  = . - — — J £ C^fo,  ' ) + 

4-  v \ + 

it4-2«+’  V 1{i  J 
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